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Abstract 

We study the dynamical critical behavior of multigrid Monte Carlo for the two 
dimensional Sine Gordon model on lattices up to 128 X 128. Using piecewise 
constant interpolation, we perform a W-cycle (7 = 2). We examine whether 
one can reduce critical slowing down caused by decreasing acceptance rates on 
large blocks by doing more work on coarser lattices. To this end, we choose a 
higher cycle with 7 = 4. The results clearly demonstrate that critical slowing 
down is not reduced in either case. 
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Typeset using REVTgX 



Computer simulations with local algorithms in statistical mechanics and lattice gauge 
theory close to a critial point suffer from critical slowing down (CSD). Besides other nonlocal 
methods, multigrid Monte Carlo algorithms were introduced to overcome this problem [1]. 

An attempt to understand why multigrid Monte Carlo is successful in beating CSD for 
some models while it does not work as well for others, was made in [2]: An analytic calcu- 
lation and analysis of acceptance rates for nonlocal Metropolis updating was performed. It 
was argued that for a critical model with Hamiltonian 7i(<f>) } CSD will occur if the expansion 
of (7i((f> + ip)) in terms of the shift ip contains a relevant (mass) term. If such a term is 
present, Metropolis step sizes s(Lb) on block lattices with increasing block size Lb have to 
be scaled down like e(Lb) ~ L^ 1 in order to obtain block size independent acceptance rates. 
This strong decrease of step sizes on large blocks was found in several models for smooth 
and piecewise constant interpolation. One of these models is the Sine Gordon model in two 
dimensions. 

The purpose of this note ist twofold: First, we want to check the prediction [2] that 
a W-cycle (cycle control parameter 7 = 2) with piecewise constant interpolation will not 
eliminate CSD in the rough (massless) phase of the Sine Gordon model. (This is a method 
that eliminates CSD in the Gaussian model.) Secondly, we want to ask the question whether 
one can circumvent slowing down caused by too small steps on large blocks by accumulating 
many of these steps randomly. If this would be possible, the acceptance problem could be 
solved by doing more work on coarser lattices. 

A constant accumulated step size on all length scales can be achieved in the following 
way: For step sizes scaling down like s(Lb) ~ Lg 1 an d a coarsening by a factor of two, 
the Metropolis step size on a next coarser grid is too small by a factor of two. Assuming a 
random- walk like accumulation of the steps, one can expect to compensate for this decrease 
by increasing the number of updates on the next coarser grid by a factor of four. This can 
be achieved by a higher cycle with cycle control parameter 7 = 4. The rule for higher cycles 
is that from an intermediate block lattice one proceeds 7 times to the next coarser lattice 
before going back to the next finer lattice. In this way 7 times more updates on each coarser 
lattice are performed. 

For a recursive multigrid algorithm, the computational effort is ~ L d for 7 < 2 d and 
~ L d logL for 7 = 2 d in d dimensions [3]. Therefore a higher cycle with 7 = 4 is practical 
for d > 2 and borderline practical for d = 2. 

The 2-d Sine Gordon model is defined on an L X L lattice A via the partition function 

Z = I J] #,exp(-^)) , (1) 

^GAo 

with the Hamiltonian 

W) = yb £ - ^) 2 - C E cos 2 ^ • ( 2 ) 

The first sum runs over all bonds in the lattice. From the point of view of statistical 
mechanics, this system can be considered as a 2-d surface in a periodic potential. The model 
exhibits a Kosterlitz-Thouless phase transition at f3 c ((). In the limit of vanishing fugacity (, 
f3 c takes the value 2/tt = 0.6366 .... For f3 > f3 c the model is in the rough (massless) phase. 
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There, the cosine-term of the Hamiltonian is irrelevant in the renormalization group sense. 
The system is critical and has the same long distance behavior as the massless Gaussian 
model. The fluctuations of the surface are given by the surface thickness 
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(3) 



where <f> denotes the average of the field over the lattice. In the rough phase, a 2 scales with 
l0gjL[4] \ 

Our simulations are organized as follows: In order to allow for high cycle control pa- 
rameters 7, we use a recursive multigrid algorithm, piecewise constant interpolation and a 
staggered coarsening with a factor of two as described in [3]. As pre-smoothing and post- 
smoothing operation, we choose a sweep of single hit Metropolis updates. The maximum 
Metropolis step size s(Lb) is scaled down like L^ 1 . Then, acceptance rates of approximately 
50% are observed on all block lattices, in accordance with the theoretical analysis of [2]. 

We study the dynamical critical behavior of the different versions of the algorithm in the 
rough phase, where the correlation length is infinite and the physical length scale is set by 
the linear size of the lattice L. Thus, we expect the autocorrelation time r to diverge with 
the dynamical critical exponent z like r ~ L z . 

The simulations were performed at f3 = 1.0, ( = 0.5. This is deep in the rough phase. 
Note that in the limit ( — > oo which corresponds to the discrete Gaussian model [5], the 
critical coupling is f3 c = 0.7524(8) [6]. Starting from an ordered configuration, measurements 
were taken after equilibration at each visit of the finest lattice. From the autocorrelation 
function for the observable A 

P4t) = {A2) _ {A)2 , (4) 

we computed the integrated autocorrelation times 

oo 

T tnt ,A = Tj + J^PAt) ( 5 ) 
Z t=l 

and the corresponding errors by a window method [7] with a self-consistent truncation 
window of 4r 8 ' nt for the energy E = L~ 2 Y,{x,y){4>x — 4> y ) 2 and the surface thickness a 2 . We 
checked that the autocorrelation functions for a 2 showed an exponential decay. 

The numerical results are given in Table I for 7 = 2 and in Table II for 7 = 4. The 
autocorrelation times Ti nt are measured in the number of visits on the finest lattice. Note 
that our runs are longer than 10 000 Ti nt (longer than 4 000 Ti nt on the 128 2 lattice). Fig.l 
shows the dependence of r 8ntj(7 2 on L for 7 = 2 and 7 = 4. For comparison, we plotted lines 
which correspond to z = 2. 

If we fit our data for the autocorrelation time of the surface thickness in the range 
32 < L < 128 with the Ansatz r mtj(7 2 = cL\ we obtain z = 1.86(4) with xV dof = °-22 for 
the W-cycle (7 = 2), and z = 1.86(4) with x 2 /dof = 4.6 for the higher cycle (7 = 4). The 
uncertainty in z is dominated by the relative error of Ti nt on the largest lattice, which is 
about 6% in both cases. We therefore estimate 

z = 1.9(1) for the W-cycle (7 = 2), 
z = 1.9(1) for the higher cycle (7 = 4) . 
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Thus, as already predicted in [2], CSD in the rough phase of the Sine Gordon model is not 
reduced by a W-cycle with piecewise constant interpolation. According to the acceptance 
analysis, we would expect the same result with smooth interpolation. Moreover, the results 
clearly show that compensating for decreasing acceptance rates on large blocks by choosing 
a higher cycle with 7 = 4 does not improve the dynamical critical behavior of the algorithm. 
We conclude that a random-walk like argumentation as stated above is not correct in the 
case of the Sine Gordon model. 
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TABLES 



TABLE I. Numerical results 


for the W-cycle (7 = 


= 2) m the 2-d bine 


Gordon model 


on£xi 


lattices in the rough phase, (3 = 


1.0, ( = 0.5. 








L statistics 


discarded 


E 


T int,E 


a 2 




4 25 000 


2 000 


0.934(4) 


0.90(3) 


0.268(1) 


0.96(3) 


8 50 000 


2 000 


0.986(1) 


0.97(2) 


0.3809(9) 


1.35(3) 


16 100 000 


2 000 


0.9956(5) 


1.04(2) 


0.4896(7) 


2.70(8) 


32 300 000 


2 000 


0.9987(2) 


1.03(1) 


0.5996(7) 


8.54(19) 


64 500 000 


2 000 


0.99945(5) 


1.04(1) 


0.7105(10) 


30.5(1.0) 


128 500 000 


4 000 


0.99966(3) 


1.04(1) 


0.8218(19) 113.7(6.9) 


TABLE II. Numerical result! 


3 for the higher cycle 


with 7 = 4 in the 2-d Sine Gordon 


model on 


L X L lattices in 


the rough phase, (3 = 1.0, ( = 0.5. 








L statistics 


discarded 


E 


T int,E 


a 2 




4 25 000 


2 000 


0.940(4) 


0.89(3) 


0.268(1) 


0.91(3) 


8 25 000 


2 000 


0.986(2) 


0.94(3) 


0.380(1) 


1.14(4) 


16 25 000 


2 000 


0.9965(10) 


0.94(3) 


0.488(1) 


1.67(6) 


32 100 000 


2 000 


0.9985(3) 


0.95(1) 


0.5997(9) 


4.15(11) 


64 300 000 


2 000 


0.99945(7) 


0.96(1) 


0.7113(9) 


14.2(4) 


128 300 000 


2 000 


0.99962(4) 


0.95(1) 


0.8213(18) 


58.2(3.3) 
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FIGURES 




FIG. 1. Dependence of the integrated autocorrelation time for the surface thickness a 2 on the 
lattice size L in the rough phase of the 2-d Sine Gordon model, (3 = 1.0, ( = 0.5. Errors are smaller 
than the symbols used. The lines correspond to z = 2. 
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